J production in two-photon collisions at next-to-leading order 
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In this paper, we report on the calculation of the cross section of J/^ plus jet inclusive production in direct 77 
collisions at next-to-leading order within the factorization formalism of nonrelativistic quantum chromodynamics 
(NRQCD). Theoretical predictions for the future e^e~ Unear collider TESLA are also presented. 



1. INTRODUCTION 

Heavy quarkonium played an important role 
in establishing the asymptotic freedom of QCD, 
as its mass is much larger than the lovif-energy 
scale of the strong interactions Aqcd- On the 
other hand, the compensation of colour in heavy 
quarkonium gives new insights into the con- 
finement mechanism, a nonperturbative process 
which is still not fully understood. 
The factorization approach based on the non- 
relativistic QCD (NRQCD) Pg provides a rig- 
oros and systematic framework for calculating 
perturbative and relativistic corrections to heavy 
quarkonium production and annihilation rates. 
This formalism implies a separation of short- 
distance coefficients, which can be expressed as 
perturbative expansions in the strong-coupling 
constant as, from the long-distance matrix ele- 
ments (MEs), which must be extracted from the 
experiment. The MEs are predicted to scale with 
a definite power of the heavy quark relative ve- 
locity V. In this way, the theoretical calculations 
are organized as double expansions in as and v. 
Within the NRQCD formalism, the colour-octet 
(CO) processes contribute to the production and 
decay rates at some level. Their quantitative sig- 
nificance was assessed for inclusive charmonium 
hadroproduction observed at the Fermilab Teva- 
tron [3]. However, it is necessary to test the uni- 
versality of the MEs in other kinds of production 
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processes, such as 77 collisions. This process was 
studied at LEP2, where the photons originated 
from hard initial-state bremsstrahlung. Hope- 
fully, the future e^e^ linear collider will be built, 
for which an additional source of hard photons 
would be provided by beamstrahlung. 
However, the leading-order theoretical predic- 
tions to be compared with the experimental data 
suffer from considerable uncertainties, mostly 
from the dependences on the renormalization and 
factorization scales and from the lack of informa- 
tion on the nonperturbative MEs. In this paper, 
we report on the calculation of the next-to-leading 
order (NLO) corrections to the inclusive produc- 
tion of J/^t mesons in 77 collisions 4 . This is 
a first step in a comprehensive programme con- 
cerning the study of charmonium production at 
NLO within the NRQCD framework. 

2. NLO CORRECTIONS TO J/* 
DIRECT PHOTOPRODUCTION 

In the following, we take into consideration the 
process 77 — + J/^P + j + X, where X denotes the 
hadronic remnant possibly including a second jet. 
We restrict the analysis to j and X free of charm 
quarks and to finite values for the J/^* transverse 
momentum pT- Two-photon processes can be 
modelled by assuming that each photon either in- 
teracts directly (direct photoproduction) or fluc- 
tuates into hadronic components (resolved pho- 
toproduction). Thus, the above mentioned pro- 
cess receives contributions from the direct, single- 
resolved, and double-resolved channels. Further- 
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more, we concentrate on the NLO calculation of 
the inclusive cross section for direct photoproduc- 
tion. 

At LO in as , there is only one partonic subprocess 
that scales like u^, namely 

ll^cc['sf^]+g. (1) 

The representative Feynman diagrams together 
with the analytic expression of the LO cross sec- 
tion for the process can be found in Ref. 
In the calculation of the production rates beyond 
LO in tts, ultraviolet (UV), infrared (IR), and 
Coulomb divergences arise and need to be regu- 
larized. One of the most convenient method to 
handle the UV and IR singularities which are 
present in the short-distance coefficients as well 
as in the MEs is dimensional regularization. We 
introduce a 't Hooft mass ^ and a factorization 
mass M as unphysical scales, and formally distin- 
guish between UV and IR poles. The Coulomb 
singularities are regularized by assigning a small 
relative velocity between the c and c quarks. 
For the extraction of the short-distance coef- 
ficients we apply the projector formalism of 
Refs. PQ. 

In the following, we study virtual and real ra- 
diative corrections separately, as their calculation 
requires different approaches. 

2.1. Virtual Corrections 

Feynman diagrams that generate the virtual 
corrections to the process can be collected into 
two classes. The first one is obtained by attach- 
ing one virtual gluon in all possible ways to the 
tree-level diagrams. They include self-energy, tri- 
angle, box, and pentagon diagrams. Loop inser- 
tions in external gluon or c-quark lines are taken 
into consideration in the respective wave- function 
renormalization constants. The self-energy and 
triangle diagrams are in general UV divergent; 
the triangle, box and pentagon diagrams are in 
general IR divergent. The pentagon diagrams 
comprising only abelian gluon vertices also con- 
tain Coulomb singularities which cancel out simi- 
lar poles in the radiative corrections to the op- 
erator / O" 



For the analytical treat- 
ment of the abelian five-point functions, we refer 
to Ref. 0. 
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Figure 1. Light-quark box diagrams contributing 
to the virtual corrections to 77 — s- cc['^S'^^-'] -I- g. 



The diagrams of the second class comprise light- 
quark loops. The individual contributions aris- 
ing from the triangle diagrams are equal to zero 
according to Furry's theorem [Hj. The box dia- 
grams, selectively shown in Fig. ^ contain UV 
and IR singularities, but their sum is finite. 
The UV divergences comprised in self-energy and 
triangle diagrams are cancelled upon renormal- 
ization of the QCD gauge couphng gs = ^/Amxl, 
the charm-quark mass m and field and the 
gluon field A. We adopt the on-mass-shell (OS) 
scheme to renormalize m, ^, and while for 
gs we employed the modified minimal-subtraction 
(MS) scheme. 

2.2. Operator Renormalization 

For a consistent NLO analysis, one should 
also consider higher-order corrections in as to 



within 



the four-quark operator [O^ ^sf^ 

NRQCD. This bare d-dimensional operator has 
mass dimension d — 1 while its renormalized ver- 
sion, which is extracted from the experimental 
data, has mass dimension 3. We thus introduce 
the 't Hooft mass scale of NRQCD, A, to compen- 
sate for the difference between bare and renor- 
malized operators. The corresponding tree-level 
and one-loop diagrams are depicted in Fig. 5 of 
Ref. 7j. Using dimensional regularization, the 
NRQCD Feynman rules in the quarkonium rest 
frame, and firstly expanding the integrands as 
Taylor series in 1/m, and afterwards performing 
the integration over the loop momentum, we ob- 
tain the unrenormalized one-loop result 
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where the subscript labels the tree-level quan- 
tity and /Lt is the 't Hooft mass scale of QCD. 
The presence of UV divergences indicates that 



the (O 



^S^^'^ ^ operator needs renormaliza- 

tion. For this, we choose the MS scheme so that 
the l/ejjv pole in Eq. ||2J) is cancelled out by the 
l/e(7V + ln{4:Tr) — 7b pole in the coefhcient of the 
counterterm ME. The term proportional with 1/v 
represents the Coulomb singularity, which com- 
pensates a similar term in the virtual corrections. 
After the renormalization of the NRQCD ME and 
cancellation of the Coulomb singularity, an IR 
counterterm at 0{as) is generated from Eq. 
which is indispensable to render the overall NLO 
result finite. This feature will be discussed in 
some detail in the next section. 

2.3. Real Corrections 

The real corrections to the process ^ arise 
from the partonic subprocesses 

7(fci) +7(^2) ^ cc[n]{p)+g{k3)+gih), (3) 
where n = ^pf) ^ 15(1) , ^ 3^(8) ^ 3pj8) ^ ^nd 
7(^1) +7(^2) ^ cc[n]ip)+q{h) + q{ki), (4) 

where n = ^S'g^'' , ^S*!*' , ^-Pj*"* . For the respective 
Feynman diagrams we refer to Fig. 2 and Fig. 3 
of Ref. 131. Note that the colour-singlet states 
n = ^S[^^ in process Q and n — ^P^' , ^S'q^^ in 
process Q are forbidden by the Furry's theorem 
[HI and colour conservation, respectively. 

Integrating the squared MEs of the processes 
iPJ and (01 over the three-particle phase space 
while keeping the value of px finite, we encounter 
IR singularities, which can be of the soft and/or 
collinear type. In order to systematically ex- 
tract these singularities, it is useful to slice the 
phase space by introducing infinitesimal dimcn- 
sionless cut-off parameters Si and Sf, which arc 



connected with the initial and final states, respec- 
tively (liij. In the case of process (jSJ we distin- 
guish soft, final-state collinear and hard regions of 
the phase space. In case of process Q we differ- 
entiate initial-state collinear and hard regions of 
the phase space. While the individual contribu- 
tions depend on the cut-off parameters 5i or df, 
their sum must be independent of them. We used 
the numerical verification of the cut-off indepen- 
dence as a check for our calculation. 

The contributions originating from process ^ 
with n ~ contain IR (soft and final-state 

collinear) singularities which cancel out singu- 
larities of the same type comprised in the vir- 
tual corrections, as required by the Kinoshita- 
Lee-Nauenberg theorem ^T]. The subprocesses 
with n = •^Py-', ■^P^*-' generates contributions con- 
taining only soft divergences. These are can- 
celled out by the IR poles leftover in the renor- 
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operator, which organise 



themselves as coefficients of the 
operators. 

According to the mass factorization theorem jl2| . 
the form of the collinear singularities present in 
the squared amplitudes of processes (0} , which are 
related with an incoming-photon leg, is universal 
and the pole can be absorbed into the bare PDF 
of the antiquark q inside the resolved photon. As 
a consequence, the real MEs acquire an explicit 
dependence on the mass factorization scale M. 
In turn, the resolved-photon contribution is eval- 
uated with the renormalized photon PDF, which 
are also M dependent. In the sum of these two 
contributions, the M dependence cancels up to 
terms beyond NLO. 

3. PHENOMENOLOGICAL STUDY 

In the following, we study the phenomenologi- 
cal significance of the process 
77 J/* + j + X at TESLA in its e+e" 
mode, with ^/s = 500 GeV and the photon 
sources bremsstrahlung and beamstrahlung co- 
herently superposed. For the energy spectrum of 
the photons we refer to Eq. (27) of Ref. ^2] and 
Eq. (2.14) of Ref. |T1], and for the designed ex- 
perimental parameters to Ref. JS] and Ref. 
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Figure 2. Differential cross section cPa/dpT dy of 
e+e- e+e-J/tp + X a.t TESLA with = 
500 GeV (a) for y = as a function of pT and (b) 
for Pt = 5 GeV as a function y. 



respectively. 

We use m = 1.5 GeV, a = 1/137.036, and the 
two- loop formula for a^s^'\^) 17 with Uf = 3 
active quark flavours and Aq^j-, — 299 MeV 
[IS]. As for photon PDFs, we employ the NLO 
set from Gluiick, Reya, and Schienbein (GRS) 
|18j . which are implemented in the fixed fixed- 
flavour- number scheme, with n/ = 3. Our de- 
fault choice of renormalization and factorization 
scales is /i = M = tot and A = to. As for J/4', 
Xcj, and V' MEs, we adopt the LO sets deter- 
mined in Ref. ^Hl using the LO set of proton 
PDFs from Martin, Roberts, Stirling, and Thorne 
(MRST98LO) 20 . 

In Fig. 12 we study d^a/dprdy (a) ioi y — 
as a function of px and (b) for p^ — 5 GeV as 
a function of y. The solid and dashed lines rep- 
resent NLO and LO results, respectively, for the 
direct photoproduction. The dotted lines corre- 
spond to LO result of the single-resolved photo- 
production evaluated with the LO versions of the 
and the photon PDFs. Notice that we 
do not consider px values smaller than 2 GeV, 
where additional IR singularities occur. 
From Fig. I3a), we observe that, the NLO re- 
sult of direct photoproduction falls off with in- 
creasing Pt more slowly than LO one. This fea- 
ture may be accounted for by the fragmentation- 
prone partonic subprocesses that contribute at 
NLO result, while they are absent at LO. Such 
subprocesses contain a gluon with a small virtu- 
ality that splits into a cc pair in the Fock state 
n = '^s[^^ . They generally provide dominant con- 
tributions at Pt 2> 2m due to the presence of a 
large gluon propagator. In single-resolved photo- 
production, a fragmentation-prone partonic sub- 
process already contributes at LO. This explains 
why the solid and dotted curves in Fig. Efa) run 
parallel in the upper pT range. In the lower pT 
range, the fragmentation-prone partonic subpro- 
cesses do not matter, as the gluon propagator be- 
comes finite and the Fock state n — '^s[^^ is al- 
ready present at LO. 

As one can notice from Fig. Hfb), at px = 
5 GeV, single-resolved photoproduction is still 
overwhelming. The two pronounced maxima of 
the LO single-resolved result may be explained 
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by the occurrence of a virtual gluon in a i chan- 
nel that can be almost coUinear with the incoming 
quark q or by the one of a virtual quark in the 
u channel that can become almost collinear with 
the incoming photon. The NLO result increases 
towards forward and backward directions due to 
the finite reminders of the initial-state collinear 
singularities that were absorbed into the photon 
PDFs. 

4. CONCLUSIONS 

The experimental verification of the NRQCD 
factorization hypothesis is of great importance 
especially for charmonium, because the charm 
quark mass might not be large enough to justify 
the nonrelativistic approximation. 
In this paper, we studied at NLO the inclusive 
production of prompt J/^P mesons with finite val- 
ues oi pt- This is the first time that an inclusive 
2 — > 2 production process was treated at NLO 
within the NRQCD factorization formalism. A 
complete NLO calculation of prompt J/'f pro- 
duction in 77 collisions allows a NLO treatment 
of photo and hadroproduction as well. This will 
provide a solid basis for the attempt to accommo- 
date the NRQCD theoretical predictions with the 
experimental data to be taken at high-energy col- 
liders, such as the Tevatron (Runll) and HERA 
(HERA-II). 
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